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ABSTRACT 

We derive one-point functions of the loop operators of Hermitian matrix-chain 
models at finite in terms of differential operators acting on the partition func- 
tions. The differential operators are completely determined by recursion relations 
from the Schwinger-Dyson equations. Interesting observation is that these gen- 
erating operators of the one-point functions satisfy VTi+oo-like algebra. Also, we 
obtain constraint equations on the partition functions in terms of the differential 
operators. These constraint equations on the partition functions define the sym- 
metries of the matrix models at off-critical point before taking the double scaling 
limit. 
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1. Introduction 

Recently much progress has been made on the matrix model formulation of 
the 2D gravity to study the non-perturbative effects^ ' and an interesting connec- 
tion with the integrable systems has been made in the double scaling limit, in 
which the size of the matrix N becomes infinite and the matrix potentials have 
critical forms.'^' In this limit, the non-perturbative results can be obtained from 
non-linear integrable differential equations, such as KdV equations. Furthermore, 

[2 — 3] 

the correlation functions satisfy their hierarchical equations. 

It has been noticed recently that the integrability of the matrix models is 
maintained even at off-critical points (finite N) before taking the double scaling 
limit. At finite N, the Lax pair, zero-curvature conditions, and infinite number of 
conserved quantities of the matrix model have been derived and related to inte- 
grable systems in more clear and direct way. The underlying integrable systems 
have been identified with ID Toda hierarchy for one-matrix model'*' which be- 
comes the KdV hierarchy in the scaling limit,'^' 2D Toda hierarchy,'*'^' and 2D Toda 

[7] 

multi-component hierarchy for the two-matrix model and for the general multi- 
matrix models, respectively. The partition functions are the 'r-functions' of these 
integrable systems. 

Next object of interest is the correlation functions of local operators. For the 
operators appearing in the action, the correlation functions are simply given by 
the derivatives of the partition functions with respect to the coefficients of the 
operators in the action. It requires, however, non-trivial analysis for the operators 
which do not appear in the action. In this paper, we derive one-point functions for 
the general local operators which are the 'loops' in the matrix models in terms of 
differential operators acting on the partition functions. These operators satisfy the 
recursion relations coming from the Schwinger-Dyson equations. We notice that 
the commutation relations of these differential operators are similar to those of the 
iyi_i_oo algebra and become exact in the continuum limit. 

One important related problem is the symmetry structure of the matrix models. 
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The Virasoro and Wp+i algebras have been conjectured for the p muhi-matrix 
models as constraint equations on the partition functions in the double scaling 
limit!* The derivation of these symmetries, however, has not been made except 
for the one-matrix model (the Virasoro algebra) and a special two- matrix model 
(the Ws algebra). This derivation may be possible if one consider the constraint 
equations of the matrix models at finite N first. Indeed, it is at finite N that the 
Virasoro algebras have been derived for the one-matrix model and for the 
multi-matrix model. In this paper, we derive most general constraint equations 
for the multi-matrix models in terms of the generators of the VFi+oo-like algebra. 
These constraint equations seem to be consistent with the conjectures made in the 
double scaling limit. 

2. Two- Matrix Model 

The partition function of the Hermian two-matrix model is given by 

Z[{tk};{sk},c]^ J VUVVe-^, 

00 q (1) 

S = Vi(C/) + V2(y) - cUV, Vi{U) = V2{V) ^J2^kVK 

k=l k=l 

Note the difference in the two potentials Vi and V2; Vi is arbitrary polynomial 

potential and V2 is with fixed order. We want to express correlation functions 
in terms of t^'s and their derivatives acting on the partition functions. These 
differential operators depend explicitly on the another parameters, s/j's. 

The most interesting loop operators in the two-matrix models are Tr(y"^C/"*). 
The one-point functions of these loops are 

Tr(v''U"'^^ = jvUVVe-^ Tr(v''U'^^ . (2) 



Prom the Schwinger-Dyson (SD) equations. 



J2 J -^uvv 



ynjjm) ^-S 



0, (X = C/,y), (m,n>0), (3) 
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one can derive two recursion relations as follows: 



m— 1 

c(Tr( = (Tr(v''U'^V[{U)j^ - J] (^Tr (^l^'^t/'^-^-i) (4) 

n-1 

c^Tr(^\/"?7^+i)^ = ^Tr(^l^"V^(y)C/"^)^ - ^^Try*Tr(^y"-*-iC/"*) ^ (5) 



s=0 



Differential operators generating the one-point functions, defined by 

W^^^h[{tk};{sk},c] ^ -c-(Tr{v-U-^)), (6) 
satisfy the recursion relation from Eq.(4), 



r=0 r=l (7) 

wS^ = j—, ■7^ = -N, for m,n>0, 

dtm Oto 

where the symbol ° is defined by {A°B)Z — A{BZ). This recursion relation can 
be rewritten in a simple form 

r^oo ^ ^^r, if r < 



where ^ ' ' ' x denotes the normal ordering. The operator J^'s satisfy C/(l) current 
algebra [Jm, Jn] = m5m+n,o {J{z) = Xlm JruZ'"^^^ = )■ Eq.(8) defines recur- 
sively the generating differential operators of the one-point functions. If there is 
no upper limit in the summation range {m + n ^ oo), it is obvious that 

w^^\z) = j2 wi['h-'^+^ = : {d,<i>y: . (9) 

m 

These infinite number of currents W^'^\z) {n — 1,2, • • •) generate the H^i+oo al- 

[13] 

gebra. For the finite values of m + n, however, the commution relations are not 



4 



exactly same as those of the VTi+oo algebra. This W^i+oo-hke algebra generates the 
one-point functions of the loops. It is remarkable that in the continuum limit the 
loop operators are given by the operators like Tr(f/*^) with 'lattice spacing' a — > 
and 'lattice size' M — > oo while keeping aM finite. Therefore, the H^i+oo algebra 
generates the one-point functions in the double scaling limit. 

For exphcit examples and later use, we write exphcit expressions for Wm , Wm , 



r=0 r>0 



m+1 m—r /m+1 m+r 

9rdsdm-r-s + XI ( X] ^sdm+r-s + ^ ^sdm+r-s ) (10) 

oo 



+ \-rtstd I irn + 2){m+l) 

' / ^ ' '^r'^'^s'^m+r+s ~r 2 
r,s=0 



where dm — d/dtm and Wm can be identified with the Virasoro generator as 
it satisfies the classical Virasoro algebra [L^, L„] = (m — n)Lm+n- 

Now consider the constraint equations on the partition function. One can 
express Eq.(5) with the one-point generating operators as follows: 

W!n^^h[{tk}; {sfc}, c] = 0, for n > 0, m > -n, 

= pFir+') - T ^w^'+:^ - T wK'^ow^i+'^ ^^^^ 

yym y^m ^"m-k Z-^^^m+k -k 

k=l k=0 

Not all of these constraints are independent. In fact, we can prove that only Wm^^s 
are independent by showing the following relation from Eqs.(4) and (11) 

m+n—l 

W!^^'^ - E Jr^^^-r: (n>0,m> -n). (12) 

r=— 00 

If Wm^Z — 0, Wm'^^^Z — are automatically satisfied. Therefore, the constraint 
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equations for the two-matrix models become 



d 

dtm 



k=l ^ 



m—k 



(13) 



As one can see in Eq.(13), the constraints are hnear combinations of the generators 
of the l^i+oo-like algebra with the coefficients of the second potential V2. In the 
continuum limit, the constraint equations are given by the algebra. Further- 

more, for a special potential V2, the operators VFm^'s may generate a subalgebra 

of the Wi+00, say, the Wn algebra. 

One can realize the VFi+oo-like algebra as a symmetry of the matrix model in 
the context of quantum field theory. In the ordinary quantum field theory, the 
symmetry can be found as infinitesimal changes of the quantum fields which leave 
the action invariant. For the matrix model, this can be done by the following 
generators Am,n, Bm,n defined by 



N r, r 

dUij [\ Ji,j 



N 



-S 



= E 



d 



(14) 



The generators A^ n and B^ n satisfy the following closed commutation relations: 



[^m,0, -Ba;,z] = kBm+k-l,l, [Am,n, Bq^i] = -n-A^^n+l-l, 
[Am,l, Bi^i] = -Bm,/+1 - ^m+l,l- 



(15) 



Note that (^Am,n[e~ ]j = (^Bm,n[e~ ]j = from the SD equations (3), which 
become the constraint equations as shown above. This realization makes it simple 
to prove the statement that only the n = constraints are independent. This 
comes from the fact that ( BmnW~^] / = can be obtained from ( Bmo\e~^] ) = 



by using the commutation relations (15). In fact, it is not difficult to see that 
{Aq^O) ^2,1, Bi Q, Bq 2} are enough to generate the constraints. 
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3. Multi-Matrix Models 



The multi-matrix models with p Hermitian matrix variables Ua have the par- 
tition function 

^ p ( p p-i 1 

Z[{tk}\ {Sa,k}. {Ca}\ = / n T^Uae-^, ^ = ^ ^ E ^«(^«) " E ^^UaUa+l \ ■ 

0=1 1 0=1 0=1 J 



We will choose the matrix potentials 



Vi(C/i) ^Y.^kVl. Va(C/a) = ^S«,fcC/,^ 



(16) 



(17) 



considering t^'s as variables and Sa^^'s as fixed parameters. The loop operators 
in the multi-matrix models are given by Tr [f/p'' ■ ■ ■ U'^'U'^'^ . Again, we want to 
express one-point functions of these operators in terms of the linear differential 
operators. 



We start with the SD equations: 



d 



N 

E / m 

N 

E / m 



d 



0, 



. . . . . . [/-.^) . . = 0, (2 < a < p - 1) 



^ ^ ij 



0. 



(18) 



Eq.(18) can be used to derive recursion relations for the one-point functions of 
the loop operators: 



ni~l 

E 

r=0 



Tr(t/;^---[/f [/fiV((C/i))), (19) 
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Tr (C/;^ • • • KltV'MKr,' • • • t/r) ), 2 < a < p - 1 (20) 



r=0 

+ (iv(v;(c/,)c/;^c/;r/ • • -c/r)). (21) 

Eq.(19) can be rewritten in the form of the recursion relations in terms of the 
differential operators 

m+n— 1 

r= — CO (22) 

We want to show that any one-point function can be expressed in terms of the 
variables t^s and their derivatives. Since Eq.(22) for ns = • • • = = with the 
identification c = ci is exactly same as Eq.(8), Wm^^^(0, • • • , 0) = Wm~^^^ of the 
two-matrix model. 

To derive other one-point functions, we consider other recursion formulae com- 
ing from Eq.(20), (2 < a < p — 1) 

yVm^^\ns, ■ ■ •,na-i,0,na+i + l,---,np) 

= -^^yVm^^\n3, Ua-l + 1, 0, Ua+l, ■ ■ ■ ,np) 

Ca-2 (23) 

+ XI -^^^m^^H'^S: • • • , ria-i, k-1, Ua+l, • • • , rip). 
k=l 

Assuming that we can express all operators W^r?''(0, ■ ■ ■ , 0, na+i, ■ ■ ■ , n^) in 
terms of t^s, we can findrf(0,---,0 ) ^a+i) ■ ■ ■ ) i^p) by repeatedly using Eq.(23). 
Therefore, we showed inductively that all the one-point functions can be found 
as differential operators of t^s acting on the partition function. Finally, if one 
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finds all the operators in terms of t^'s, one can find the constraint equations on 
the partition function from Eq.(21). Among others, the case of = gives the 
following equations: 



(rig, • • • , + 1, 0)Z[{tk}; {sa,k}, {ca}] = 0, 
Wff (ns, ■■■,np-i + 1, 0) = Wff (na, • • • , n^-i + 1, 0) 

- E ^>^S^ (ns, • • • , ^ - 1), (p > 4) 

fe=i Vi 



and we must treat carefully the index n in Eqs.(23) and (24) for p < 3. 
We apply above general analysis to the three-matrix model. Defining 



they satisfy the following recursion relations: 



m+n—l 



>v«(/ + 1) = -ciww ,(0 + ^wff (0. 

k=2 



(24) 



y\^t-l\l)Z[{tk}] = -cM(Tv(iyV"C/-)), (25) 



(26) 



As explained above, from 'Wrn\^) = one can find W,'n''(l)'s from the second 

equation and win*''(/)'s from the first one. Continuing this step, one can find all 
Wm'*(0'^- Finally, the constraint equations come from Eq.(24). For an explicit 
example, consider the following potentials: V2{V) — V2V'^ + vsV^ and V3(VF) = 
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W2W'^ + W2,W^\ From Eq.(26), one can find the explicit expressions: 



Cl c| 

m r) oo 

(1) = E m;:^inUi^ + 5]rt,w« ,(1), 

r=0 r=l 

m+1 oo 

r=0 ^ r=l 

Wff(2) = -ciWW ,(1) + + 

Cl C-|^ 



(27) 



and the constraint equations are 



»fl,-^>vi:'(i)-5£ipw(;)(2) 



^2 ^2 



0, (2^ 



where are given in Eq.(lO). 

4. Discussions 

In this paper, we computed one-point functions of the multi-matrix models in 
terms of the differential operators acting on the partition functions. The operators 
arc completely determined by the recursion formulae, derived from the SD equa- 
tions and generate the 14^i+oo-likc algebra. Furthermore, we derived the constraint 
equations on the partition functions using these operators. Since the partition 
functions are the 'r functions' of the 2D Toda hierarchies,'*'*^'^' this means the one- 
point functions as well as the partition functions are completely determined by the 
integrable systems and symmetry structures. Our method can be generalized to 
the multi-point functions. Again, the generating operators are determined by the 
recursion relations which comes from the SD equations. 

It is also possible to consider the s^^'s as variables such that one can introduce 
another differential operators for the two- matrix model. For the potentials Vi{U) — 
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EktkU\V2{V) = Ek^kV\ one can define 



wt^l^Z[{tk},{sk}] ^ -c-(iv(c/-y-)), 



(29) 



where both Wm and W^^ satisfy the recursion relation Eq.(8). The constraint 
equations are just 



Another interesting point we want to mention is that the constraint equations 
for the multi-matrix models like Eq.(28) have very similar form as those of the two- 
matrix models Eq.(24). The coefficients of the second potential of the two- matrix 
model can be decided by those of the multi- matrix potentials. The correspondence, 
however, is not quite exact. There appear some terms in the constraint equations 
of the multi-matrix models which do not exist in those of the two-matrix model. 
This observation reminds us of the recent conjectures that all the multi-critical 
points can be achieved by the two-matrix model.'^*' If the extra terms at finite 
N are suppressed in the double scaling limit, our observation can be a proof of 
this claim. Related to this and other motivations, it would be very interesting 
to consider the double scaling limit of our formalism. Our discovery that the 
correlation functions are generated by the H^i+oo algebra acting on the 'r-functions' 
of 2D Toda hierarchy seems to be consistent with the results in the double scaling 
limit in that the correlation functions of the one-matrix model are given by KdV- 
hierarchy equations and that the r-function of the p reduced KP-hierarchy satisfy 
the Wi+p constraint equations. 

Recently, there have been several papers which mention the Wi+oo algebra. 
Our Wi+oo algebra is different from that of rcf. [15] in that the latter comes from the 
higher order terms under the change of M to M -|- 5M. Therefore, this constraints 
exist even for the one-matrix model. Our VFi+oo constraints exist only for the 



m—n 



(n+1) 




Z[{tk},{sk}] = 0- 



(30) 
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multi-matrix models and will have direct connection with the Wn algebra structures 
conjectured in [8] in the double scaling limit. Also, the VTi+oo algebra appears from 
the KP hierarchy in the double scaling limit This is a direct p ^ oo limit of the 
Wi+p constraint equations considered in [9]. It would be interesting to consider 
the p ^ oo limit of our result to understand these results from the matrix model 
point of view. 

Note: While typing this paper, we received a paper [17] from Y.-X. Cheng where 
constraint equations for the two- matrix model like Eq.(30) have been obtained. 
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